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This paper considers the steady state queueing properties of pedestrians who are delayed by 
traffic on a main road. It is assumed that the size of the gap between two successive cars determines 
the probability that an individual pedestrian or a group of pedestrians will cross. The technique of 
imbedded Markov chains yields many new results as well as some previously found, by combina- 
torial methods by Tanner in 1951. 

1. Introduction 

There has been a recent increase of interest in the problem of the delay to traffic at a stop 
sign. This problem was first suggested by data analyzed by Adams in 1936, [l]. 1 Adams derived 
an expression for the mean delay to a pedestrian caused by a stream of traffic on a main highway, 
in which the successive gaps (in time) are independent random variables, each with a negative 
exponential distribution. It was assumed that each pedestrian waits for a gap which is larger than 
a prescribed safe interval before attempting to cross. Garwood, in 1941 [2], took up the problem 
in a modified form. He dealt with the problem of driver delay at an intersection in which a traffic 
signal is partially operated by the main stream traffic and partially by the crossing traffic. When a 
crossing car arrives at the intersection it actuates the signal so that the light will change to green 
either at a fixed time M after its arrival, or at some earlier time if the main stream is sufficiently 
light. Garwood's analysis assumed a negative exponential distribution for headways. Garwood 
derived an exact distribution for the delay. The next discussion of the delay problem was given 
by Raff, [3], who also assumed a negative exponential gap distribution and a gap acceptance func- 
tion given by a step function, i.e., cross if the gap is greater than a given T c . The same problem 
was taken up in much greater detail at about the same time by Tanner [4] who derived a large num- 
ber of results for the statistics of the size of the queue of pedestrians. He assumed that the 
pedestrians' arrival distribution was Poisson, as well as that of the traffic, and calculated such quan- 
tities as the expected number of pedestrians crossing between successive cars, the size of the queue 
at a random time, and others which will be reproduced in the present paper. The basic assump- 
tion made by Tanner is that of a step function distribution for the gap acceptance function both for 
an individual pedestrian and for the group. We might note that the process described by Tanner, 
and indeed that which is described below in the following manuscript, is a bulk service process in 
which a stationary state is always reached for all gap distributions observed in moderate traffic. A 
later contribution to the pedestrian delay problem was made by Mayne, [5], who showed, by fairly 
involved calculations, how to calculate the delay to a single pedestrian when the gap distribution 
was other than Poisson, but otherwise assuming a step gap acceptance function. 

All of the work mentioned above was carried out by means of combinatorial arguments, which, 
except in the case of a negative exponential gap distribution, were quite complicated. Recent 
work on delay problems has proceeded from the observation that the passage of a car on the main 
highway is a regeneration point, and the techniques of renewal theory can be used to solve delay 
problems in a way which almost completely eliminates combinatorial reasoning from the solution. 
As an example of papers which make use of renewal-theoretic techniques we cite [6-9]. Because 
the solution of delay problems has now been freed of the complications of the earlier analyses it is 
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natural to attack more realistic situations. Thus, in the first paper on renewal theory as applied to 
the delay problem, not only was it possible to relax the restriction of a negative exponential gap 
distribution (this had already been done by Mayne) but it was also possible to take into account a 
gap acceptance function which was not a step function. Specifically, one can define a probability 
a(t), such that an individual will cross the highway with probability a(t) when the gap which con- 
fronts him is of duration t. Actual measurements of a(t), both on pedestrian and vehicular traffic, 
indicate that a{t) is a continuous function of t rather than a step function, [10-12]. A typical form 
found in measurements on the vehicular gap acceptance probability is 

a(t) = 0^t^T c (11) 

= l-exp[-p(t-T c )] t^T c . 

The theory which has been developed and partially checked against experiment applies to a single 
driver at a stop sign confronted by a single lane of homogeneous traffic. Very little is known about 
driver behavior when confronted by a multilane highway, although it is not difficult to extend the 
theory to the case where the waiting driver considers the situation in each lane independently. 
One can also solve for the moments of the delay distribution for a simple model of nonhomogeneous 
traffic in which there are TV types of cars and the gap distribution between successive cars depends 
on the types involved, [6]. 

While the single car, or single pedestrian delay distribution is a handy function for the traffic 
engineer to know, it has little bearing on the queueing properties of the system. This is a difficult 
problem to deal with when the queues are considered to be those of cars, who must necessarily 
form into a line, and may be released in groups whose size depends on the size of the gap in traffic 
on the main highway. In this situation, where it is assumed that the second car in line takes a non- 
negligible time to reach the head of the line when the first car has crossed, the traffic on the main 
highway may indeed be so heavy as to cause an unlimited buildup of a queue on the feeder road. 
In the pedestrian queue first studied by Tanner all of the pedestrians in the queue see the gap simul- 
taneously and can cross simultaneously. It is this feature that insures that except for pathologic 
and obviously unrealistic gap distributions, the queue will have steady state properties. 

It is the purpose of this paper to reconsider the pedestrian queueing problem, eliminating the 
restriction to a negative exponential gap distribution and using a gap acceptance probability which 
is not necessarily a step function. It will be seen that by using the techniques associated with 
imbedded Markov chains (first introduced into this problem area by Oliver and Bisbee, [7]), the 
resulting calculations become simple, because of the particular form found for the transition 
matrix, when it is assumed that the gap acceptance function for a group of pedestrians is the same 
as that for a single pedestrian. We shall calculate the probability generating functions and the 
first two moments for the number of pedestrians in queue at each regeneration point, the number 
in queue at a random time, and the number crossing after a random vehicle. 

2. General Theory 

We begin by considering the following specification of the situation to be studied. We 
assume: 

1. Poisson arrival of pedestrians, with arrival rate X. 

2. General stationary headway distribution in which successive gaps are independent random 
variables described by a density (p(t). 

3. The gap acceptance function for each pedestrian is a(t). 

The first and simplest problem for pedestrian queueing is the study of the imbedded Markov 
chain, in which we take as the regeneration points for the process, the times of departure of cars 
on the highway from a point directly in front of the group of pedestrians. We shall calculate the 
elements of the transition matrix P = (p W m) in several situations. The element p nm represents the 
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probability that if there are n pedestrians in queue at one regeneration point, there will be m 
pedestrians in queue at the succeeding regeneration point. To calculate p nm we must be slightly 
more precise about our definition of regeneration points. It will be assumed that if a car passes in 
front of the queue at time t, any pedestrians who choose to cross will do so at a time £ + 0. 

The problem which corresponds to those treated by previous authors is that in which either a 
pedestrian crosses on arrival at the intersection (i.e., he is not delayed at all) or else the whole 
group of pedestrians in queue decides to cross in a body when a new gap presents itself. The gap 
acceptance function for the group is taken to be that for the individual. We shall now calculate 
the pnm for this situation. 

There are two cases to be examined. When n is greater than m the transition n-+m can take 
place only as n— »0^ra, that is, first the queue empties entirely and then m or more pedestrians 
arrive of which m choose to join the queue. When m is greater than or equal to n the transition 
n-^m can take place in two ways. The first is n^O— >m as already discussed, and the second is 
n—*m. In this second possibility the pedestrians in queue choose not to cross the road as the car 
passes and m — n pedestrians arrive who choose to join the queue. 

Since we are assuming Poisson arrival of pedestrians, the arrival times can be considered to 
be uniformly distributed in the interval (0, t) for any t, [14]. Thus, the probability that a pedes- 
trian will see an initial gap of (x, x + dx) units of time will be dx/t. Therefore, the probability of 
immediate crossing by an arbitrary arrival in this interval is a*(t), where 



'Jo' 



a*(t) = - a(x)dx, (2.1) 



a*(t) is the average of a(t) taken with a uniform weight over the interval. The matrix element 
Pnm for n > m is therefore 

= j °°a(t)e-W-«*W [kt{l ~^* (t))]m <p(t)dt, n>m. (2.2) 



The derivation of this relation is as follows: The transition n—>0 (the entire queue crosses with 
the passage of a car) takes place with probability a(t); k+m pedestrians arrive in the gap of t 
units, with probability e~ kt {kt) k+m l(k-\- m)l; of these m remain in queue, where m is a random variable 
with a binomial distribution with probability 1 — a*(t). Notice that the second line can be sum- 
marized by assuming a Poisson arrival rate on a reduced time scale denned by 



T(t) 



= t[l -«*(*)]= J (l-a(u))du. (2.3) 

Jo 



When m is greater than or equal to /i, another term is needed in addition to that given in eq (2.2). 
This term is needed to cover the case of the direct transition n~* m and is 

^oJo (k + m — n)\ \ m — n J 

f* (\T(t)) m ~ n 

= [1 - ««] V , t e-™<p(t)dt (2.4) 

Jo (m — »)! 
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where T(t) is defined in eq (2.3). The transition matrix P is therefore 

Ao + e Ai + d A 2 + € 2 A 3 + e 3 . . . 

A Ai + e A 2 + ei A 3 + e 2 . . . 

A Ai A 2 + € A 3 + ei . . . 

A Ai A 2 A 3 + e . . . 



(2.5) 



/here 1 



Jo 



k\ 



e k = j [l-a(t)]<p(t)^^e' kT dL 



(2.6) 



The gross properties of this Markov chain are easily deduced. Let us denote by 6(n) the vector 
of state probabilities after the passage of the nth car, where 



The vector 6{n) satisfies 



or 



fKn) = Mn) 9 0i(n),. . .)• 

9(n+l) = 0(n)P 
d(n) = 0(O)P n . 



(2.7) 

(2.8) 
(2.9) 



(2.10) 



Because of the simple form of P shown in eq (2.5), eq (2.8) is most conveniently written in component 
form k 

fe(n + l) = A* + y0Xra)€fe-j 

which suggests the use of generating functions. Defining the generating functions 

Vn(s) = %6 k (n)s k , e(s) = %e k s k 



A'=0 



A-=0 



we find 

and by iteration that 

It is readily verified that 



A(5)=£a***, 

k=0 



V n + 1 (s) = A(s) + e(s)Vn(s) 

1 - € n {s) 



y n (s) = 6 n (s)y (s) + A(s) 



l-e(s) 



Ms) 

e(s) 



= [ " a(t)<p(t)e- kT « ni - s) dt 

Jo 

= J X [1 - a(t)Mt)e- mt)(1 - s) dt. 

Jo 



(2.11) 

(2.12) 
(2.13) 

(2.14) 



2 Although the derivation of p j = Aj + €j is not quite correct, since one cannot speak of a group of persons choosing to cross the road, the result is nonetheless 
correct, since Aj + €j is independent of a(t). 
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It would be difficult to solve explicitly for all of the components of 0(n). However, it is not 
difficult to find an expression for the moments. Let us denote by fik(n) the kth ordinary moment, 
that is 

fjL k (n) = ym k m (n). (2.15) 

»1=0 

Noting that 

e(l) + A(l)=l (2.16) 

we find for the first two moments 

/xi(n) = mi(°°) + &0) IMO) ~ Mil 00 )] 

Mn) = lLx(n) + €»(1) [/x 2 (0) - ja,(0)] + \Z^$ j " (Wf<p(t)dt 

l-e»(D 



+ 2 



H,(™) + ne n -\l)(fii(0)-fJLxM) 



L 



kT(l - a(t))<p(t)dt (2.17) 



where 
and 



L 1-ed) 

e(\)={ X [l-a(t)] tf >(t)dt (2.18) 

J() 

I *KT{t)<p(t)dt I °° \T(t)<p(t)dt 

M»)=— — -==7= < 2 - 19 ) 

l-e(l) a{t)<p(t)dt 

Jo 



The equilibrium value of the second moment, /M 00 ) is given by 






M«(~) = /Xl(O0)+ l- e (l) + Q 



(2.20) 



By these methods we can determine the distribution of the several parameters discussed by 
Tanner and by Mayne. The probability that there are n pedestrians in queue after a random 
vehicle is 0n(°°). Hence the expected number of pedestrians crossing after a random vehicle 
passes is 

Yn0nM\ a(t)<p(t)dt = f i l (cc) J a(t)<p(t)dt. (2.21) 

n=0 Jo Jo 

The kth moment of the distribution is similarly equal to /x/c(°°) I a{t)(p(t)dt. 

Jo 

The probability that there are no pedestrians in queue when the nth vehicle passes is 

6 (n) = ^ w (0) = e*(0)^o(0) + A(0) \Z^ (2.22) 



where 

e(0) 



= €o=] a{t)<p(t)e- xnt) dt 

= A = r[l-a(t)]<p(t)e- kTit) dt. 
Jo 

At equilibrium {n = &>) we have the result 



(2.23) 
A(0) = 



•w-i^l-i^;- 



233 



The probability density for the amount of time that any single pedestrian is delayed will be denoted 
by £l(t), and can be immediately taken from the results of reference [6], since each pedestrian can 
be regarded as not being influenced by the presence of any others. Let us denote the Laplace 
transform of any function, say F(t), by the same function with argument s and an asterisk, e.g., 

F*(s)=f C °e- gt F(t)dL (2.25) 

Jo 

Then it is shown in [6] that 

a*(s) = j a(t)<p (t)dt + j^T^j f a(t)<p(t)dt (2.26) 

where 

&t) = <p(t)[l-a(t)] 

(2.27) 
t ) o(t) = <Po(t)[l-a{t)] 

and where <po(t) is the density of the first gap given by 

<Po(t)= <p(x)dxH x<p(x)dx. 

Expressions for the first two moments of the delay density are given in [6] and are repeated here for 

reference: 

— T7oo 
t^iqoi+T- 1 t?h 

1 — T?10 

9 9 2 < 2 ' 28 ) 

Vo2 1-ijio (l-^o) 2 



where 



Voj=\ t j £ > o(t)dt=\ tJ(f (t)[l — a(t)]dt 
Jo Jo 

VU = I X t^(t)dt= I °°t?<p(t) [l-a(t)]dt. 

Jo Jo 



(2.29) 



Another quantity of interest that can be derived from the present analysis is the distribution 
of the number of pedestrians crossing between two consecutive vehicles. Let 77/,- be the probability 
that k pedestrians cross between two successive vehicles. In solving for rrk we must take two cases 
into account; either the group of pedestrians initially in queue (when the first car passes the group) 
cross together, or it does not. The probability that k or more pedestrians arrive in a time t (the 
time delay between the arrival between two successive cars) and that exactly k of these cross, is 

fe-^-J k + m ) [«W[1-*W -i^^r^^-r^MWiM^Mt (2.30) 

where T(t) is defined in eq (2.3). Hence, considering separately the two cases mentioned above, 
we have, for the probability that k pedestrians cross assuming that the group initially in queue 
crosses. 



k r °° 

■ (1) =20K°°) <* 

j=0 J 



*) E l^-^- (2.3D 



In the second case, where it is assumed that the group in queue does not cross we have, 

W 2) = f"[l -a(t)]<p(t) [Xto *j t)] * e -x**o<fr, (2.32) 
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In this situation the size of the initial queue is obviously irrelevant. Finally ttk- is just the sum 

7T k = 7T k il) -\-TT k {2 \ (2.33) 

The moments can be derived from the generating function 

rr(s)=yir k s k = a(%KO^oo(s)e- x ' a * ( ' )(1 - s Vft + [l-a(t)](p(t)e- xta *^-^dt. (2.34) 

k=o J$ Jo 

The first two moments of this distribution, denoted by ri\ and n 2 are 

r oo r oo /» oo 

/ii=jLti( 00 ) at(t)(p(t)dt+ kta*(t)(p(t)dt = \ t<p(t)dt 
Jo Jo Jo 

roc ^oo 

n 2 = M 0O )\ a(t)<p(t)dt + 2^(00) l \tQL*(t)a{t)<p{t)dt 
Jo Jo 

+ [ [\to*W]VW* + ^. (2.35) 

Similar calculations enable us to find the probability that there are m pedestrians in queue at 
a random time, not necessarily just after the passage of a car. As before we shall assume that the 
event "passage of a car" is a regeneration point for the process, and calculate a set of transition 
probabilities, {k nm }, where A. nm is the probability that there are m pedestrians in queue at a random 
time, given that there were n in queue at the last regeneration point. For the present calculation 
we assume that the gap in which the random time finds itself is of duration t, and that the last re- 
generation point occurred y units of time ago, where clearly y ^ t. 

We first calculate the probability density for the duration of the gap, t, in which the instant 
findi itself. Let w(x)dx be the probability for the occurrence of a regeneration point to occur in 
(x, x + dx). Let p(t, y)dt be the probability that a point at y is covered by a gap whose duration 
satisfies t ^ G < t + dt. We shall only be interested in the limit y^> oo. The probability that any 
gap has duration G such that t ^ G < t + dt is just <p(t)dt. However, we must account for the fact 

that the last regeneration point lies in the interval 3 y, y— t. This probability is just I w(x)dx, 

Jy-t 

hence the probability density for the size of gap including a random time point is 



p{t,y) = <p(t) 



I w(x)dx - I w{x)dx +0(1) (2.36) 



in the limit of y^ 00 . Under quite general conditions on <£>(£) which are always met in practice it 
can be shown by the methods of renewal theory, [13], that 



fy 

w(x)dx ~ 
Jo 






(2.37) 



(2.38) 



in the limit of y - ^ 00 , where 

V\ = t(f(t)dt. 
Jo 

Hence we finally find 

p(«,oo) = -^. (2.39) 

This result has been derived in a slightly different manner by Oliver and Jewell, [15]. 



3 We tacitly assume that y is chosen so large that y> t. A more complete result which does not use this assumption is easy to prove, but in what follows we are 
only interested in the equilibrium result. 
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We now define a function b nm (y,t) by 8 n m(y, t)dy= probability that there were n pedestrians 
in queue at the last regeneration point which occurred in (t — y, t — y+ dy), and that m pedestrians 
are presently in queue. The quantities which we seek, the \ n m are related to d nm {y, t) by 

knm=\ p(t)dt\ -J§nm(y,t) (2.40) 

since the time of observation is uniformly distributed over the gap duration. An expression for 
$nm(y, t) can immediately be written down by exactly following the argument which led to eq (2.5). 
The final result is 



- n > m 



8„ m (y, t) = a(t)e-^y-» [XT+iy ; t)] ' m^n 

nl 



(2.41) 



(m— n)l 



where 



KT + (y, t) = k\ (l-a(u))du. (2.42) 

Jt-y 



Hence, finally, k n m can be expressed as the sum of two terms: 

\nm=A m n>m 

=A m + B m -n m^n 
where, combining eqs (2.38), (2.40), and (2.42) 



(2.43) 



A tt = -— ( a(t)<p(t)dt f [KT+^Dfe-^y'^dy 
n\V\ J o J o 

B r = -7— I (l~a(t))<p(t)dt I [\T + (y, t)] r e~ XT+ ^ t] dy. 
r\v x Jo Jo 



(2.44) 



Let us denote by c m the probability that there are m pedestrians in queue at a random time. 
This probability is expressible in terms of the ^ n ( oc ) as 

oc m 

Cm = 2 0nMKm = A m + £ n (™)B m - n . (2.45) 

n = n = o 

Since this equation is in convolution form we can profitably write down the generating function 
for c m , call it C{s), and relate it to the generating functions of A m and B m , A{s) and B(s) respectively. 
In this way we are led to 

C(s) = A(s) + t-^J- B(s) (2.46) 

where we have used eq (2.13) for the generating function of the Oni 00 ). The generating functions 
A{s) and B(s) are explicitly 

A(s) = — | ' a(t)<p(t)dt \ e -^(y,t)(i-s) dy 
Vi J o Jo 

B(s) = - [~ {I- dt))<p{t)dt [ e- kT ^^- s) dy. (2.47) 

V\Jo Jo 
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_ A 
m~- 



-2 -.*' 

m l = m -t - 



For the first two moments we find by differentiation 

— \ <p{t)dt |V + (y,^y+ i ^ [™ t{l-(x{t))ip(i)dt 

^1 Jo Jo ^1 Jo 

'■\*<p{t)dt [ [T+(y 9 t)]*dy 

Jo Jo 

+ ^^ f X (1 - <*))tft)<ft [' r+(y, t)dy 

^i Jo Jo 

+ &M f x t{ i- a{t))ip{t)dL (2.48) 

^i Jo 
When a(£) is a step function, 

a{t) = t<T c 

= 1 t>T c (2.49) 

the double integrals can be reduced to single integrals by introducing the explicit form for T + {y, t): 

T + (y,t) = t>T c , t-y>T c 

= T c -t + y t>T c , t-y<T c 

= y t<T c . (2.50) 

When this form for T + (y, t) is substituted into eq (2.47), we obtain the results 

2viJo 2v x v\ Jo 

A 2 f T c A 2 ? 13 

m 2 = m + —\ tWt)dt + ^<P(T c ) (2.51) 



V\ J o 



^1 



fhere 



<&(t)= [ <£>(x)d*. (2.52) 

3. Some Applications 

We shall now apply some of the preceding general results to the specific distributions 

{t) = ae -*t^ m = 0,l,2,... 
^ ml 

a(t) = t<T c 

= l-e-W~V t^T c . (3.1) 

When /3 tends to infinity this form of a(t) tends to the step function 

a(t) = t<T c 

= 1 t^T c (3.2) 
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which has been used extensively in previous treatments of this subject. It is easily verified that 
for this choice of a(t) we have for a*(t): 

a*(t)= t<T c 

, (3.3) 

= j t [J3t-pT c -l + e-w-V} t>T c 

and 

T(t) = t[l-a*(t)] 

= t t<T c 



(3.4) 



= h\ + pT c -e-w-V] t^T c 



A quantity which frequently occurs in the calculations is e(l), which for the present set of models is 

e(l) =]<p m (iil -a(t)]dt 

Jo 

J ml Jt c n\ 

= 1 - G m (aT c ) + (— ^ j eFcGJLio- + p)T e ] 



(3.5) 
where 

GJM = er*\l + x+^+. . .+—X (3.6) 



When /3 = <x >, the expression for e(l) becomes 

e(l) = l-G m (oT c ). (3.7) 

An expression for the quantity 

[\T{t)<p m {t)dt 

which appears in the formula for ^ti(°°) (of eq (2.19)) is 

f \T{t)<p m {t)dt = (m+ 1)- [1 - G m + i((tT c )] 
Jo cr 

+ 1 (1 + fST c )G m (<rT c ) - 1 ePc (^g) m + 1 G m [(<r + ffl t \ (3.8) 

When (3 = !X this expression reduces to 

f mt)Vm{t)dt={m+ 1)- [1 - G m + i(o-T c )] + \T c G m (crT c ). (3.9) 

Jo O" 

As a specific example of the use of these formulas, let us specialize to the case m=0, i.e., 

<p(t) = ae-^ (3.10) 

where a is the expected number of cars passing per unit time. Then it is easily found that 

e(D=l-^^e-^o (3.11) 

and 

^ k $ + j)h-p&) (3 - i2) 
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Figure 1. Plot of normalized expected queue length as a function of 
normalized traffic intensity. 



is the equilibrium value of the expected number of people in queue at a regeneration point. When 
|8 = oo we obtain the value given by Tanner for the case of a negative exponential headway 
distribution: 

iLtl (oo)=A (e crT c _ 1) (313 ) 



We find in general for /3 = °°, the value of ^i(°°) is 

> 
/xi(oo) = xr c + (/n+l) ' 



CT 



1 — G m + i(oT r ) 



Gm(crTc) 



(3.14) 



Several curves showing /xi(oo)/(Xr c ), the normalized expected number of pedestrians in queue at 
the passage of a random car are shown in figure 1. As cr— > oo, i.e., the mean gap between vehicles 
decreases to zero, we find 



Mi(«)=xr c [^p^^c+o(i)]- 



When (r-^Owe find 



Mi(») ~ kT e 1 + 



[' 



(<TT c ) m ^ 
m\(m + 2) 



(3.15) 



(3.16) 



The equilibrium probability that there are no pedestrians in queue when a random car passes 
is given in eq (2.24). When <p{t) is given by (p m (t) of eq (3.1) and a(t) is a step function #o(°°) is 



698-888 0-63-4 
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Figure 2. Plot of probability of an empty queue as a function of 
normalized traffic intensity. 



0o(°°) = 



e-VcGnkrTc) 



/ (j. \m+l 
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In particular when m = this expression reduces to that derived by Tanner, [4], 



0o(«>) = 



ke {a + x)T c + cr 



(3.17) 



(3.18) 



Figure 2 contains plots of 0o(°°) for \T C = 1 and different values of m. When a — > o° eq (3.17) yields 
the approximate result 
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(3.19) 



and when o - — »0, or equivalently, traffic is light, eq (3.17) reduces to 

1 , G m (\T c ) 



o (oc) = € -kt c ^ _ ((rTc)m + 1 j_i 



+ i)! (xr c ) m+i (xr c ) 



si: 



+o((o-r c r +i ). 



(3.20) 



In the first approximation we see that the probability that there are no pedestrians in queue does 
not depend on the vehicular traffic, but only on the pedestrian arrival rate. 

The expected number of pedestrians crossing between two cars can be calculated from eq 
(2.35) When <p(t) is of the form of eq (3.1), then detailed calculation leads to the result 



7ii = (m+ 1) — • 
cr 



(3.21) 
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The expected size of the queue at a random time is given in the first line of eq (2.47). When 
(fk{t) and a(t) are of the form given in eq (3.1) we find 



=4 { —^ [1 - C* +2 (oT c )] +^\zG k {aT c . 



2 a- 



2(k+l) 



+ [l-G fc+ ,(crr c )] 



Xr c + (A:+1) : 



l-G k+ i(<rT c ) ] 
■ G k (o-T c ) J 



When cr— »co this expression becomes approximately 

__ (£+l)!Ae°- r c 
m ~ *{aT c T 

and for cr- -»0 the expression goes into 



+ 0(1) 



(3.22) 



(3.23) 



m = kT c [l + 0(aT c )]. 
For the Poisson distribution we set k = and obtain the exact result 



as derived by Tanner. 



m=^(e- T c-l-(TT c ) 



(3.24) 



(3.25) 



4. Comments and Generalizations 

One of the weakest points of the preceding analysis is the assumption that a group of pedes- 
trians will have the same gap acceptance probability as a single pedestrian. If we assume that 
there are, instead, an ensemble of gap acceptance probabilities {a n (t)} then the form of P is not the 
same as that shown in eq (2.5), and it is necessary to discuss the resulting Markov chain ab initio. 
In particular it is no longer possible to have as simple a relation between generating functions as 
that given in eq (2.12). The more complete form for P is easily shown to be 



where 
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e v = f "(1 - ait))e-^l - a*(t)] lm ~" ? ® V " <p(t)dt. 



(4.2) 



Attempts to solve the steady state equations corresponding to the transition matrix of eq (4.1) have 
been unsuccessful. It is no longer clear in this situation that the queue will reach a steady state, 
although this will evidently be the case when a\(t) ^ az{t) ^a 3 (t) ^ . . . since we have seen that a 
steady state is reached when all of the a(t) are equal, and the inequality provides for faster servicing 
of the arriving pedestrians. 
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There is one case of differential gap acceptance probabilities that can be solved. Let us 
assume that 

a l (t) = P(t),a 2 (t) = a 3 (t)=. . .=a(t) (4.3) 

that is, that group behavior will differ from individual behavior, although the gap acceptance 
probability for a group of greater than one will not depend on the size of the group. We shall 
examine the steady state equations only, i.e., the solution to the equation 0(°°)= 0(°°) P. When these 
equations are written out in component form (omitting the infinity in the argument) we find 

n 
On = A n + ^ Qr€n-r + d\y n - (4.4) 

r=0 

where Aj and e, are just the functions defined in eq (2.6), and 



y n =\ *>(t)e- x TO) -a(*)] 



\\T) n (KT) n - 1 
n\ (w-1)!. 



dt (4.5) 



70= ( <p(t)e-* T \fi(t)-a(t)]dt 
Jo 

and 



T(t)=f (1-P(u))du. (4.6) 

Jo 



The equation satisfied by the generating function is 

oc x 

where "^(5)— V 6 n s n and y{s)= Vy/. The right side of this equation still contains the unknown 

value of 0i, but this can be eliminated by differentiating both sides and setting s equal to zero. 
With this procedure we find 

l = A,(l- eo ) + Aoe 1 

Wl (l-6 ) 2 — *(1 -€„)->€, (4 ' 8) 

which can then be substituted back into eq (4.7) to give the full expression for the generating 
function. 

When it is assumed that the first r of the ctj(t) differ, but the remaining aj(t) are all equal, then 
r of the 0's appear in the expression for the generating function. These can be found by inverting 
an rxr matrix. We will not pursue this line of inquiry any further since no experimental data exists 
which is relevant to this matter. 

Another shortcoming of the model, even in the somewhat more general formulation of the last 
paragraph, is the fact that either none or all of the pedestrians in queue are assumed to cross at 
the same time. One can formally remedy this deficiency by defining a matrix {a nr (t)) such that 
a nr {t) is the probability that r out of n pedestrians will accept a gap of duration t. However, de- 
tailed calculations show that the resultant transition matrix is such that the stationary problem 
cannot be solved by generating function methods, and so far the solution to the problem is unknown. 

It is possible by our methods to treat the queueing problems with bunched arrivals, where the 
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bunches arrive in accordance with a Poisson process, but the number in each bunch is a random 
variable with a known distribution. Again, detailed calculation for this situation is hardly infor- 
mative without any supporting data. 



I thank Alan Goldman for his careful reading of this manuscript and for several suggestions 
which have been incorporated into the final draft. 
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